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: ABSTRACT 

^ I We have measured the present accretion rate of roughly 800 low-mass (~ 

^1 1 — 1.4 Mq) pre-Main Sequence stars in the field of Supernova 1987A in the 

' Large Magellanic Cloud (LMC, Z ~ 0.3 Zq). It is the first time that this fun- 

' damental parameter for star formation is determined for low-mass stars outside 

' our Galaxy. The Balmer continuum emission used to derive the accretion rate 

■ - - ' positively correlates with the Ha excess. Both these phenomena are believed 

to originate from accretion from a circumstellar disk so that their simultaneous 
detection provides an important confirmation of the pre-Main Sequence nature 
of the Ha and UV excess objects, which are likely to be the LMC equivalent of 
Galactic Classical T Tauri stars. The stars with statistically significant excesses 
are measured to have accretion rates larger than ~ 1.5 x 1O~^M0 yr~^ at an 
age of 12-16 Myrs. For comparison, the time scale for disk dissipation observed 
in the Galaxy is of the order of 6 Myrs. Moreover, the oldest Classical T Tauri 
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star known in the Milky Way (TW Hydrae, with 10 Myrs of age) has a measured 
accretion rate of only 5 x 10"^'^ M^lyr^ i.e. 30 times less than what we measure 
for stars at a comparable age in the LMC Our findings indicate that metallicity 
plays a major role in regulating the formation of low-mass stars. 

Subject headings: stars: formation, pre-main-sequence, galaxies: Magellanic 
Clouds 



1. Introduction 

The processes at play during star formation determine much of the appearance of the 
visible Universe. The shape of the stellar Initial Mass Function (IMF) and its normaliza- 
tion (the star-formation rate) are, together with stellar evolution theory, key ingredients in 
determining the chemical evolution of a galaxy and its stellar content. Yet, our theoretical 
understanding of the processes that lead from diffuse molecular clouds to stars is still very 
tentative, as many complex physical phenomena concur in producing the final results. While 
clear variations in the star-formation rate are observed in different regions of the Milky Way 
and in external galaxies, with their histories showing bursts and lulls (e.g., Tolstoy 2000), 
the observational evidence for (or against) variations in the IMF is often contradictory (see 
the review of Scalo 1998, or the Gilmore vs Eisenhauer debate in "Starbursts: Near and Far', 
2001). Yet, variations in the IMF can dramatically alter the chemical evolution of a galaxy 
(e.g., Wyse 1998). 

From an observational standpoint, most of the effort has traditionally been devoted to 
nearby Galactic star-forming regions, such as the Taurus-Auriga complex, Orion, etc. If 
this, on the one hand, permits one to observe very faint stars at the best possible angular 
resolution, on the other it is achieved at the expense of probing only a very limited set of 
initial conditions for star formation (all these clouds have essentially solar metallicity, e.g., 
Padget 1996). 

Studying the effects of a lower metallicity on star formation is also essential to under- 
stand the evolution of both our own Galaxy, in which a large fraction of stars were formed 
at metaUicities below solar, and what is observed at high redshifts. As a matter of fact, the 
global star formation rate appears to have been much more vigorous (a factor of 10 or so) 
at 2; ~ 1.5 than it is today (Madau et al 1996, and subsequent incarnations of the so-called 
"Madau plot"). At that epoch the mean metallicity of the interstellar gas was similar to 
that of the Large Magellanic Cloud at present (LMC, e.g. Pei et al 1999). This fact makes 
the study of star forming regions in the LMC especially important for the understanding of 
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galaxy evolution. 

With a distance modulus of 18.57±0.05 (see the discussion in Romaniello at al 2000), the 
LMC is our closest galactic companion after the Sagittarius dwarf galaxy. At this distance 
one arcminute corresponds to about 15 pc and, thus, one pointing with a typical imaging 
instrument comfortably covers almost any star forming region in the LMC (10 pc see, e.g., 
Hodge 1988). In particular, the field of view of about 2.7' x 2.7' of the WFPC2 on board 
the HST corresponds to 37 pc x 37 pc, and leads to the detection of several thousands of 
stars per pointing (e.g., Romaniello et al 2002). The LMC is especially suited for stellar 
populations studies for two additional reasons. First, the depth of the LMC along the line of 
sight is negligible, at least in the central parts we consider (van der Marel and Cioni 2001). 
All of the stars can, then, effectively be considered at the same distance, thus eliminating 
a possible spurious scatter in the Color-Magnitude Diagrams. Second, the extinction in its 
direction due to dust in our Galaxy is low, about E(B — V) ~ 0.05 (Bessell 1991; Schwering 
1991) and, hence, our view is not severely obstructed. 

There is currently a widespread agreement that low mass stars form by accretion of ma- 
terial until their final masses are reached (e.g. Bonnell et al 2001, and references therein). As 
a consequence, the accretion rate is arguably the single most important parameter governing 
the process of low-mass star formation and its final results, including the stellar Initial Mass 
Function. Cround and HST-bascd studies show that there may be significant differences 
between star formation processes in the LMC and in the Calaxy. For example, Lamcrs et al 
(1999) and dc Wit et al (2002) have identified by means of ground-based observations high- 
mass pre-Main Sequence stars (Herbig AeBe stars) with luminosities systematically higher 
than observed in our Galaxy, and located well above the "birthline" of Palla and Stabler 
(1991). They attribute this finding either to a shorter accretion timescale in the LMC or 
to its smaller dust-to-gas ratio. Whether such differences in the physical conditions under 
which stars form will generally lead to differences at the low mass end is an open question, 
but Panagia et al (2000) offer tantahzing evidence of a higher accretion also for LMC low 
mass stars. 

In this paper we present the first measurement of the accretion rate onto low-mass pre- 
Main Sequence stars outside of our Galaxy. The data and the reduction process are presented 
in the next section, while the detection of Balmer continuum excess and its conversion into 
an accretion rate are described in section 3. Section 4 is devoted to deriving the stellar 
parameters of our sample, i. e. the stars' masses and ages. Finally, the conclusions are drawn 
in section 5. 
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2. Observation and data reduction 

The field of SN 1987A in tlie LMC was repeatedly imaged over the years with the 
WFPC2 on-board the HST to monitor the evohition of its Supernova remnant. We have 
taken advantage of this wealth of data and selected from the HST archive a uniform dataset 
providing broad-band coverage from the ultraviolet to the near infrared, as well as imaging in 
the Ha hne. The log of the observations we have used is reported in Table 5. A description of 
the camera and its filter set can be found in Biretta et al (2001) . All images are centered with 
the Planetary Camera chip on SN 1987A (asooo = 05 : 35 : 28.26,^2000 = -69 : 16 : 13.0), 
but have different position angles on the sky, resulting in complete coverage of an almost 
circular region of 130" (about 30 pc) in radius. 

The data were processed through the standard Post Observation Data Processing Sys- 
tem pipeline for bias removal and fiat fielding. In all cases cosmic ray events were removed 
combining the available images after accurate registration and alignment. 

The plate scale is 0.045 and 0.099 arcsec/pixel in the Planetary Camera and in the 
three Wide Field chips, respectively. We performed aperture photometry following the pre- 
scriptions by Gilmozzi (1990) as refined by Romaniello (1998), i.e. measuring the fiux in a 
circular aperture of 2 pixels radius and the sky background value in an annulus of internal 
radius 3 pixels and width 2 pixels. Due to the undersampling of the WFPC2 Point Spread 
Function, this prescription leads to a smaller dispersion in the Color-Magnitude Diagrams, 
i.e. better photometry, than PSF fitting for non-jittered observations of marginally crowded 
fields (Cool and King 1995; Romaniello 1998). Photometry for the saturated stars was recov- 
ered by either fitting the unsaturated wings of the PSF for stars with no saturation outside 
the central 2 pixel radius, or by following the method developed by Gilliland (1994) for the 
heavily saturated ones. The flux calibration was done using the internal calibration of the 
WFPC2 (Whitmore 1995), which is typically accurate to within 5% at optical wavelengths. 
The spectrum of Vega is used to set the photometric zeropoints ( VEGAMAG system) . 

In the following, as a measure of the overall photometric accuracy, we will use the mean 
error in five broad bands defined as: 

^ _ ^ ^(mF336w)^ + (^(mF439w)^ + ^(mpsSSw)^ + ^(wiF675w)^ + ^(mpsHw)^ ^^^^ 
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2.1. From colors to luminosity and temperature 

Once the observed fluxes of tlie stars are carefully measured, we derive their intrinsic 
properties, i.e. luminosity and temperature, as well as the extinction caused by the interven- 
ing interstellar dust along the line of sight, using the prescriptions developed by Romaniello 
et al (2002). The intrinsic stellar parameters and their associated errors are derived with a 
minimum technique by comparing the observed magnitudes to the ones expected based 
on the theoretical stellar atmosphere models of Bessell et al (1998) computed in the HST- 
WFPC2 bands using the IRAF synphot task. In order to cope with the effects of interstellar 
dust we have used the the extinction law appropriate for this region of the LMC (Scuderi 
et al 1996). Let us stress here that, by convolving the theoretical spectra with the filter 
sensitivity curves provided in IRAF, we ensure that the observations are faithfully modeled. 
In particular, the red leak that affects the WFPC2 F336W (U-band-like) filter is properly 
taken into account. 

The dereddening method is extensively described in Romaniello et al (2002), but let us 
briefly summarize it here: 



stars for which both E{B — V) and Tg// can be simultaneously derived are selected 
according to their photometric error (65 < 0.1) and their location in the Qubi vs {U — I) 
plane, where Qubi is a reddening- free color defined as: 

Q,,,^(U-B)-^^(B-I) 

This color-based selection is aimed at solving the possible non-monotonicity of broad- 
band colors with temperature (see, for example, Allen 1973, p207). The stars for which 
E{B — V) and Tg// can be derived simultaneously turn out to be hotter than 10,000 K 
or between 6,750 and 8,500 K. 

Also, a star's location in the Qubi vs {U — I) plane provides a first guess of its tem- 
perature and reddening; 

for the stars selected in step 1, E{B — V), Tgff, L and their associated uncertainties 
are derived starting from the first guesses by performing a minimum multi-band fit 
of synthetic colors from Bessell et al (1998) to the observed magnitudes; 

for each star for which E{B — V) and Tgff cannot be derived simultaneously because of 
its intrinsic temperature and/or too large errors, the reddening is set as the mean of the 
ones of its 4 closest neighbors with direct reddening determination. The corresponding 
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rms is used as an estimate of the uncertainty on the adopted value of E{B — V). 
The effective temperature, luminosity and associated errors are, then, derived from a 
minimum multi-band fit. 

In the case of the field discussed here, there is, on average, one star with direct E{B—V) 
determination every 13 arcsec^. Of course, it is possible that a few stars have, in reality, 
extinction values significantly different from the local mean, but the global effect is 
negligible. 

The errors on E{B — V) and Tg// are computed from the maps and propagated to 
the luminosity. A full discussion on the errors is reported in Romaniello ct al (2002), but it 
is important to keep in mind here that the procedure outlined above does not introduce any 
systematic errors on the derived stellar parameters. 

Panagia et al (2000), using this same dataset, have measured the age of the massive 
stars in the field to be 12 ± 2 Myrs. Also, they have identified several hundreds pre-Main 
Sequence stars through their Ha emission. We will now measure, for the first time outside 
our Galaxy, the accretion rate on these low mass pre-Main Sequence stars. 

3. Meeisuring the accretion rate 

The idea that the strong excess emission observed in some Galactic low-mass, pre-Main 
Sequence stars (T Tauri stars) is produced by accretion of material from a circumstellar disk 
dates back to the pioneering work of Lynden-Bell and Pringle (1974). The excess luminosity 
is, then, related to the mass accretion rate. In particular, the Balmer continuum radiation 
produced by the material from the disk as it hits the stellar surface has been used as an 
estimator of the mass infall activity (see, for example, GuUbring et al 1998, and references 
therein) . 

3.1. The U-band excess 

The presence of a rich population of stars with U-band excess in the field of SN 1987A is 
clearly detected in Figure 1, where we plot as gray dots the dereddened {m-p2,z6w,Q — ?7t.f439W,o) 
vs (mF439w,o ~ "^F8i4W,o)) i-^- (Uq — Bq) vs (Bq — Iq), colors for stars with good overall 
photometry (^5 < 0.1). The F336W and F439W filters bracket the Balmer jump and an 

excess in the (mF336W,o ~ ''^F439W,o) color translates in an excess emission in the Balmer 
continuum. The typical errorbar, computed as the mean of the photometric errors of the 
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stars with 0.6 < (mF439w,o — itif8UW,o) < 1-2, is shown as a cross. For reference, the 
expected locus from the stellar atmosphere models of Bessell et al (1998) for Z = 0.3 Zq and 
log{g) = 4.5 is shown as a sohd line. At the end of this section we will discuss the influence 
of chemical composition and surface gravity on our results. 

The reddening and intrinsic stellar parameters for the stars with possible U band ex- 
cess were computed excluding the F336W magnitudes from the multi-band minimum fit 
described in section 2.1. This is because the excess itself can make the observed F336W 
flux significantly different from that of a photosphere, hence yielding incorrect results for 
the stellar parameters. For example, the roughly 600 stars under the Main Sequence noted 
by Panagia et al (2000) in this same field have too high a temperature precisely because 
the U band was not excluded from the fit. As detailed in Romaniello et al (2002), however, 
it is not possible to rehably derive both Tgff and E{B — V) without the U band. For the 
stars with possible U band excess, then, we have used the reddening value from their 4 near- 
est neighbors with a direct determination of E{B — V) and derived the stellar parameters 
accordingly. 

The presence in Figure 1 of stars with too blue a (mF336w,o ~ "^f439W,o) color given 
the observed (mF439w,o — ttifsuw^) color is apparent at (mF439w,o — ttifsuw^) — 0.9. To 
quantify the U-band excess we consider the distribution of the (mF336w,o — ''Wf439W.o) color as 
shown in Figure 2. There, the solid line shows the color distribution of the stars within the 
color interval 0.6 < (^Tipdsow.o ~ '"^F814W,o) < 1-2. The FWHM expected from the measured 
photometric and dcrcddcning errors is shown by the horizontal errorbar. It nicely reproduces 
the width of the right-hand side of the observed histogram, while it severely underestimates 
the one on the left (the dashed fine is the histogram produced by the reflection of the red 
part of the distribution relative to its peak value). Hence, we take the resulting symmetric 
distribution, which has a half-power width of about 0.25 magnitudes, as an indication of 
the broadening caused by photometric and dereddening errors. This is justified by noticing 
that the peak of the observed (mF336W,o — ?71f439W,o) distribution occurs where the models 
atmospheres predict it to occur ((771f336W,o~''^F439W,o) = —0.27, thin vertical line in Figure 1), 
which provides an independent check on the dereddening procedure we have used, ruling out 
any significant systematic errors on the inferred dereddened colors. An excess of stars with 
F336W emission as compared to a normal photosphere, i.e. blue (mF336W,o — '^F439W,o) 
colors, is evident (hatched area). 

The location in the HR diagram of the stars with and without U excess is shown in 
Figure 3. Also indicated are the theoretical Zero Age Main Sequence and 12 Myr post- 
Main Sequence isochrone from Brocato and Castellani (1993) and the pre-Main Sequence 
isochrones by Siess et al (1997) for Z — 0.3 Zq, the metallicity appropriate for the young 
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population of the Large Magellanic Cloud (Hill et al 1995; Geha et al 1998). The typical 
uncertainties on the luminosity and temperature for the stars with U excess, computed as 
the mean of the uncertainties on the individual stars, is shown as a cross. If interpreted as 
pre-Main Sequence stars, then, the position of the stars with U excess in the HR diagram 
of Figure 3 is consistent at a with an single age of 12 Myr, the age inferred from the upper 
Main Sequence (see also Panagia et al 2000). We will return on the age dating of these stars 
in section 4. 

We can now take a step forward and compute the excess luminosity in the F336W filter, 
which, as we will see in section 3.3. is directly related to the mass accretion rate onto the pre- 
Main Sequence star. This is accomplished by computing the difference between the measured 
dereddened flux in the band and the one predicted by the stellar atmosphere models (Bessell 
et al 1998) for a photosphere with the effective temperature of the star. We have derived this 
latter quantity by fitting the observed BVRI magnitudes to the same model atmospheres of 
Bessell et al (1998): 

-^F336W,exc — 47r AAf336W {D^ /F336W,obs — ^1 -^F336W,mod) (2) 

where AApssew = 370 A is the width of the F336W filter (Biretta et al 2001), D = 51.8 kpc 
is the distance to this LMC field (Romaniello at al 2000), i?* the stellar radius and /F336W,obs 
and FF336W,mod ^re the observed and model flux densities, respectively. The stellar radius, 
like the effective temperature, is obtained by fitting the observed BVRI magnitudes to the 
same model atmospheres of Bessell ct al (1998). The U magnitude was excluded from the fit 
because it can be significantly contaminated by the excess, hence yielding incorrect results 
for the stellar parameters. 

The distribution of the resulting excesses in solar units is shown in Figure 4. The solid 
line represents the observed distribution, which is clearly skewed toward positive values. If 
neither excess emission nor observational errors were present, the distribution would be a 5 
function centered on 0. Therefore, we consider the broadening toward negative values as a 
measure of the combined photometric plus dereddening errors. By mirroring the histogram 
for I/F336W,o about (dashed histogram in Figure 4) one can identify the stars with excess 
emission as those populating the bins where the observed distribution exceeds the mirrored 
one (hatched region). An inspection to Figure 4 shows that for i^F336W,exc 

> 0.035 Lq the 

contamination due to excesses induced by random errors becomes negligible. 

These are, of course, the same stars with too blue a (mpssew.o ~ ^iF439W,o) color for 
the (mF439w,o — ?7iF8i4W,o) ouB highlighted in Figure 1. As shown in Figure 5, the Balmer 
continuum excess ivF336W,o displays a well defined correlation with the (mF336W,o — ?7t.f439W,o) 
color, albeit with a quite large scatter rms = 0.04 Lq. As expected, the correlation is in the 
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sense that the more negative the color, the higher the inferred F336W excess. 

The excess emission was computed assuming for all of the stars a metallicity appropriate 
for the young LMC population {Z = 0.3 Zq, e.g. Hill et al 1995; Geha et al 1998), in the 
sense that the observed F336W fluxes were compared to models of that metaUicity. This 
simplification does not have any impact on our results, as we expect the contamination from 
metal-poor [Z < 0.1 Zq, i.e. intrinsically brighter in the F336W), field subgiants to be 
of the order of only a few percent (Cole et al 2000). In addition, the stellar atmosphere 
models of Bessell et al (1998) indicate that the difference in the F336W luminosity between 
Z = 0.3 Zq and Z = 0.1 Zq is of the order of 0.01 Lq in the relevant range of effective 
temperature (5,500 K < Tg// < 8,000 K, see Figure 3) and for i?* = 1.5 Rq, the typical 
radius of these stars. This value is much smaller than the detection threshold we have set 
at -LF336w,exc ^ 0.035 Lq. The models of Bessell et al (1998) also indicate that the effect of 
surface gravity on photospheric emission in the F336W band is even smaller than the one 
due to metallicity. 

In conclusion, the 765 stars with an excess lager than about 0.035 Lq are bona fide 
objects with F336W excess when compared to a stellar photosphere. 

3.2. A sanity check: the correlation between U and Ha emission 

In addition to the U-band excess. Galactic Classical T Tauri stars are also observed 
to have Ha emission, which is also thought to be linked to the accretion process from the 
circumstellar disk (e.g. Calvet et al 2002). As such, a correlation is to be expected between 
these two quantities. 

In Figure 6 we plot the excess emission in the F336W filter versus the one in Ha: 
mF336w,obs is the observed magnitude, r«F336W,mod is the photospheric one from the models 
of Bessell et al (1998) and the (mF675w — TTi-pesm) color measures the Ha equivalent width 
(e.g. Romaniello et al 2003). The correlation between U and Ha activity is apparent. 

A high statistical correlation between the two quantities is confirmed by Spearman's 
coefficient p (see, e.g., Conover 1980). This is a non-parametric test sensitive to correlations 
between two random variables. The correlation coefficient is 7.9, which implies a probability 
of less than 10~^ that the two variables are uncorrelated. Whereas the quality of the data 
does not allow to derive the detailed shape of the correlation, the existence of the correlation 
itself is proven with a very high statistical significance. This provides a very important sanity 
check as to the pre-Main Sequence nature of these stars. 
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A potential source of concern here is chromospheric activity that would produce excess 
emission both in the Balmer continuum and in Ha, without being linked to the presence of a 
circumstellar accretion disk. As we will show in section 4, the stars we deal with in this paper 
have about 1 Mq (see Figure 8). The total chromospheric emission of a solar- type star on the 
Subgiant Branch corresponds to a luminosity of the order of 9.4 x 10~^ Lq (Ulmschneider 
1979; Pasquini et al 2000). Even if the flux were all concentrated in the spectral region 
covered by the F336W filter, its contribution would still be negligible compared to the 
threshold we have set of 3.5 x 10~^ Lq. Moreover, it is clear from Figure 6 that all the 
objects with significant F336W excess also have (mpersw— "7.f656n) ^ 0.15, i.e. EW(Hq;) > 3, 
which is an upper limit to what can be expected due to chromospheric activity (Frasca and 
Catalano 1994). 

In conclusion, then, the case under study the chromospheric emission does not contribute 
in any significant way to the observed Balmer continuum excess. 



3.3. From U-band excess to mciss accretion 

Now that we have detected and characterized the stars with a statistically significant 
excess in the U band we can compute the corresponding accretion rate. The following 
equation: 

.^.^lM^(:_M (3) 

relates the bolometric accretion luminosity to the stellar and disk parameters. Rin is the 
inner radius of the accretion disk. Its value is rather uncertain and depends upon the details 
of the coupling of the accretion flow to the magnetic field. Following GuUbring et al (1998) 
we adopt Rin — 5 R* for all of the stars. 

GuUbring et al (1998) also provide an empirical relation between the accretion luminosity 
Lace and the excess luminosity in the Johnson U passband compared to a stellar photosphere: 

log = 1.09 log + 0.98 (4) 

This, in turn, has been transformed to the WFPC2 F336W fllter by Robberto et al (2004) 
using theoretical models for the disk emission that reproduce GuUbring et al (1998) 's relation: 
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log (^^^ = 1.16 log (^I^E^^^^ + 1.24 (5) 

Combining equations (3) and (5) leads to a relation that links M, the accretion rate in 
Mq yr~^ to the excess luminosity in the F336W filter i^F336W,ea;c' 

log(M)[Me yr-] = -6.70 + 1.16 log (^^^E^^^ + log - log (6) 

3.4. The mass accretion rate 

According to equation (6) in order to derive the accretion rate M we need the following 
three quantities: 

• -^F336W,ea;c which is Computed as described in section 3.1 and whose resulting distri- 
bution is shown in Figure 4; 

• the stellar radius R^, which is calculated from the luminosity and temperature deter- 
mined for every star using the method developed by Romaniello et al (2002), i.e. by 
fitting the observed colors to the ones predicted by stellar model atmospheres, and 
adopting a distance modulus to the LMC of 18.57 (Romaniello at al 2000); 

• the stellar mass iVf* that we estimate by comparing a star's location in the HR diagram 
with evolutionary tracks computed for Z = 0.3 Zq (Siess et al 1997). 

The distribution of the derived accretion rates is shown in Figure 7 for the stars with 
-^F336W,ea;c > 0.035 Lq, i.e. thosc for which the excess can be measured with high reliability. 
As it can be seen, the pre-Main Sequence stars detected in the field of SN 1987A through 
their U-band excess have accretion rates larger than ~ 1.5 x IO^^Mq yr^^. 

The apparent low cutoff at low values of M is not real, but, rather, is a consequence of 
the threshold we have used to identify stars with a significant Balmcr continuum excess. In 
other words, there might be stars in the field with smaller accretion rates, but they would 
be lost in the measurement errors. Moreover, Galactic T Tauri stars are known to show 
variability at many wavelengths with timescales as short as a few days (see, for example, the 
review by Bertout 1989). It is clear, then, that, since we only detect those stars with the 
strongest excesses, our method tends to pick up stars when they are the most active. 
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4. Ages and masses of the pre-Main Sequence stars 

In section 3.1 we noticed that a visual inspection to the HR diagram indicates that, if 
interpreted as pre-Main Sequence objects, the stars with U excess are consistent with an age 
of roughly 12 Myrs. We can now refine that estimate by comparing each star's luminosity 
and temperature to those predicted as a function of mass and age by the pre-Main Sequence 
evolutionary models of Siess et al (1997) ior Z — 0.3 Zq. 

The resulting age distribution of the stars with ivpssew.eic > 0.035 Lq is shown in the 
left panel of Figure 8. The peak of the distribution (14 ± 2 Myr) is in very good agreement 
with the age of the massive stars (M > 6 Mq) in this same field (12 ± 2 Myrs, Panagia et 
al 2000). The spread in the histogram is consistent with being produced by photometric 
and dereddening errors only, which, for these stars, are of the order of 20% in the effective 
temperature (see also Figure 3). At this stage, then, we cannot assess if the stars with 
significant U band excess show an intrinsic spread in age or not. It is worth noticing that, if 
we take the comparison with the Siess et al (1997) isochrones at face value, almost no stars 
have ages lower that about 5 Myrs. 

The mass distribution of the same stars with significant Balmer continuum excess is 
shown in the right panel of Figure 8. The mean mass of the accreting stars is 1.25 Mq, with 
an rms dispersion of 0.2 Mq. 

The mean value of the accretion rate for the stars with significant Balmer continuum 
excess is ~ 2.5 x 10~^Mq yr~^ (see Figure 7). If we assume that this is a steady-state mean 
value, then these stars still have to accrete a significant fraction of their mass (10-20%) 
before settling to their final Main Sequence value. Of course, this is surely an overestimate 
for two main reasons. As we have already noticed, we tend to pick up the stars when their 
activity is maximum. In addition, accretion activity is observed to slow down with time in 
Galactic pre-Main Sequence stars (e.g. Calvet et al 2002), breaking the simphstic steady-state 
hypothesis. 

5. Summary and conclusions 

We have identified about 800 stars with statistically significant Balmer continuum ex- 
cess in the field of SN 1987A in the Large Magellanic Cloud. This excess positively correlates 
with HcK emission, as derived from the comparison of broad and narrow-band photometry. 
Whereas the quahty of the data does not allow to derive the detailed shape of the correlation, 
the existence of the correlation itself is proven with a very high statistical significance (Spear- 
man's test returns a probability of less than 10~^ that the two variables are uncorrelated) . 
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Both the Balmer continuum and the Ha excesses are well above the levels expected from 
chromospheric activity. These facts lead us to interpret these objects as pre- Main Sequence 
stars. 

In this framework, given their location in the HR diagram and the emission in the Balmer 
continuum and Ha line, the objects in our sample are the equivalent of Galactic Classical 
T Tauri stars. This is the first time that such a fundamental parameter as the accretion 
rate is measured for this class of objects in a galaxy other than our own Milky Way. Doing 
so in the Large Magellanic Cloud provides a unique opportunity to sample astrophysical 
conditions not represented in local star-forming regions, but that were common in the early 
Universe. 

When interpreted as pre-Main Sequence stars, the comparison of the objects' location 
in the HR diagram with theoretical evolutionary tracks allows one to derive their masses 
(~ 1 — 1.4 Mq) and ages (~ 12 — 16 Myrs). At such an age and with an accretion rate 
in excess of ~ 1.5 x 10~^Mq yr~^, these candidate pre-Main Sequence stars in the field of 
SN 1987A are both older and more active than their Galactic counterparts known to date. In 
fact, the overwhelming majority of T Tauri stars in Galactic associations seem to dissipate 
their accretion disks before reaching an age of about 6 Myrs (Haisch et al 2001; Armitage et 
al 2003). Moreover, the oldest Classical T Tauri star know in the Galaxy, TW Hydrac at an 
age of 10 Myrs, i.e. comparable to that of our sample stars, has a measured accretion rate 
some 30 times lower than the stars in the neighborhood of SN 1987A (MuzeroUe et al 2000). 

The situation is summarized in Figure 9, adapted from MuzeroUe et al (2000), where 
we compare the position in the age-M plane of the stars described in this paper with that 
of members of Galactic star-forming regions. An obvious selection bias that affects our 
census is that we only detect those stars with the largest Balmer continuum excesses, i.e. 
highest accretion rates. There might be stars in the field with smaller accretion rates, either 
intrinsically or because they were observed when the accretion activity was at a minimum, 
which fall below our detection threshold. This selection effect is rather hard to quantify, 
but it is clear that the locus of the accreting stars that we do detect in the neighborhood of 
SN 1987A is significantly displaced from the one defined by local pre-Main Sequence stars. 

There are essentially two ways to reconcile the position of the LMC point in Figure 9 
with the mean locus observed in the Galaxy: significantly reduce the age of the stars in the 
SN 1987A region and/or significantly decrease their inferred accretion rate. Neither option, 
however, seems viable. In fact, while the uncertainties on either quantity can be quite large 
for any given star, the LMC point reflects the mean values of several hundred stars, thus 
making the random errors, for all practical purpose, negligible. Let us now consider the 
systematic errors. 
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If we move the LMC point horizontally in Figure 9, in order to occupy the same area 
as Galactic stars, the LMC stars would have to be younger than about 4 Myrs or less, the 
age of the oldest Galactic star with a measured accretion rate of 1.5 x 10~^Mq, yr^^. This 
would require our estimate of the ages of the LMC stars to be systematically wrong by a 
factor of three or more. Such a large shift is rather implausible on accounts of two main 
considerations. First, Romaniello et al (2002) have shown that the dereddening technique 
we have adopted here does not introduce any significant systematic errors on the stellar 
parameter, hence, given a set of pre-Main Sequence evolutionary tracks, on the derived age. 
Second, the age of 14 Myrs inferred from the peak of the histogram in Figure 8 agrees very 
well both with that of the most massive stars in the field (12 ± 2 Myrs, Panagia et al 2000) 
and that of the progenitor of SN 1987A (10-12 Myrs see, for example, Scuderi ct al 1996, 
and references therein). These age estimates are, of course, completely independent from 
that of the low- mass stars we are discussing here, lending support to it. 

The other possibility is that, given the age of about 12-16 Myrs, we have overesti- 
mated the accretion rate onto the candidate pre-Main Sequence stars in the surroundings 
of SN 1987A. A factor of about 30 reduction is required in order to make our value for the 
LMC to agree with TW Hydrae, the oldest accreting T Tauri known in our Galaxy (see 
Figure 8). According to equation (6) this corresponds to reducing the Balmer continuum 
excess -Lpsaew.exc by a factor of almost 20 (for comparison, this would be so small as to be 
even smaller than the bin size of the histogram in Figure 4). Such a reduction, however, 
is clearly incompatible with the observed distribution of the excesses, in particular with its 
asymmetry, which marks the stars with a Balmer continuum excess beyond statistical errors 
(see also Figure 2). 

A potential source of concern is the effect on the derived accretion rate of the presence of 
circumstellar dust in the immediate surroundings of the candidate pre-Main Sequence stars 
(fiared disks, remnants of the cocoon they were formed out of, etc.). This extra source of 
reddening is not accounted for by our method and the candidate pre-Main Sequence stars 
could be affected by a higher value of E{B — V) than mean one from the 4 closest neighbors 
with a direct determination, i.e. have bluer intrinsic (mFssew.o ~''^F439W,o) colors. According 
to the data in Figure 5, however, this would correspond to increasing the inferred Balmer 
continuum excess and, consequently, the accretion rate, further confirming the difference 
with what is observed in the Galaxy. 

Finally, let us notice that the accretion rate for several of the Galactic stars in Figure 8 
was derived with the same technique we have used in this paper, i.e. the relation between 
Balmer continuum excess and accretion rate by GuUbring et al (1998). As such, they are 
affected by the same issues on the cahbration of this relation as the LMC stars presented 
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here (if any). The discrepancy between Galactic stars and those in the neighborhood of 
SN 1987A in the LMC, then, appears to be real. Of course, spectroscopic data on these stars 
would be highly desirable to better characterize their properties. 

A higher accretion rate in the LMC than in the Galaxy was also suggested for A and B 
spectral type pre-Main Sequence stars by Lamers et al (1999) and de Wit et al (2002). The 
combination of this result with ours points towards a significantly higher accretion activity 
over a factor of about 6 in mass for pre-Main Sequence stars at lower metallicity. 

Wc warmly thank Luca Pasquini for enlightening discussions on stellar chromospheric 
activity. The detailed comments and suggestions from an anonymous referee are gratefully 
acknowledged. 
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"^F336W.O ~ "^F439W.O 

Fig. 1. — Locus in the (mpsaew.o — "^f439w,o) vs (mF439w,o — '"^F814W,o) plane of the stars with 
good overall photometry, i. e. a mean error in the broad bands between F336W and F439W 
smaller than 0.1 mag (see equation 1). The theoretical locus from the models of Bessell et al 
(1998) for Z = 0.3 Zq and log(g') = 4, 5 is shown as a full line. The thin vertical line marks the 
position of the peak of the histogram in Figure 2. The horizontal dashed hnes highlight the 
region used to select stars with excess F336W emission (0.6 < (^f439W,o — ^F814W,o) ^ 1-2, 
see text). The presence of stars with excess F336W emission is evident inside this box 
at (mF439w,o — TTipsuwfi) — 0.9. The typical uncertainties for these stars is shown by the 
errorbar. 
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Fig. 2. — Distribution of (mpssew.o — ''^F439W,o) color for stars with 0.6 < (mF439w,o — 
^F8i4W,o) < 1-2 (see Figure 1). The fuU hne is the observed histogram while the dashed one 
is the reflection of its red part about the peak (thin vertical line), which also corresponds 
to the expected color of the turnover point in Figure 1. The asymmetry of the (mF336w,o — 
?7T.F439W,o) toward blue colors, i.e. F336W excess, is apparent and is also highlighted by the 
hatched area. The FWHM of the distribution expected from the measured photometric and 
dereddening errors on the (mpssew.o — "^f439W,o) color is shown as an horizontal errorbar. As 
it can be seen, it nicely reproduces the part of the histogram to the right of the peak, while 
it severely underestimates the width of one to the left. 
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log(T.„) [K] 

Fig. 3. — HR diagram displaying the position of the stars with U excess (black dots) overlaid 
on the general stellar population (gray squares) found in the WFPC2 field. Luminosities 
and temperatures for the stars with excess were computed excluding the F336W magnitude 
from the fit to the model atmospheres and adopting for each star the mean E{B — V) value 
of its 4 closet neighbors. The typical uncertainties on the luminosity and temperature for 
the stars with U excess, computed as the mean of the uncertainties on the individual stars, 
is shown as a cross. For reference, we also plot the theoretical Zero Age Main Sequence, 
with the position of stars of various masses marked on it, and a 12 Myr post-Main Sequence 
isochrone {Z = 0.3 Zq, Brocato and Castellani 1993). Also shown are 5-20 Myrs pre-Main 
Sequence isochrones (Siess et al 1997), again computed for Z — 0.3 Zq. 
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Fig. 4. — Histogram of the measured flux excess in the F336W fiher (full line). The stars 
with an excess emission statistically beyond the observational error are highlighted in the 
hatched area (see text). The contamination due to spurious emission becomes negligible 
above about 0.035 Lq and we set this as our detection threshold. 
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Fig. 5. — Balmer continuum excess vs dereddened (mpssew.o ~ "^f439W,o) color for the stars 
with good overall photometry (S^ < 0.1). The result of a linear regression is displayed as a 
solid line. As expected, the two quantities show a well defined correlation. 
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Fig. 6. — F336W vs Ha excess for our sample of candidate pre-Main Sequence stars in 
the field of SN 1987A (gray dots), mpsaew.obs is the observed magnitude, mpssew.mod is the 
photospheric one from the models of Bessell et al (1998) and the {mpQ^^^^ — mpeseN) color 
measures the Ha equivalent width (e.g. Romaniello et al 2003), as shown on the right vertical 
axis. The filled squares represent the median (mpersw — ^F656n) value in the (mF336W,mod — 
mF336w,obs) bins marked by the rectangles, whose vertical extent includes 66% of the stars 
in each bin. The high statistical correlation between the two quantities is apparent and is 
confirmed by the high value of Spearman's coefficient p (see text). Note that this threshold 
also excludes stars with normal chromospheric activity (3 A or (mpeysw ~ "^f656n) > 0.15, 
Frasca and Catalano 1994). 
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Fig. 7. — Accretion rates as derived from the excess flux in the F336W filter (see equation 6) 
for the stars for which the excess can be measured with certainty (i^Fsaew.exc > 0.035 Lq, 
see Figure 4) . 
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Fig. 8. — Age {left panel) and mass [right panel) distributions of the stars with highly 
rehable U-band excess (ivF336W,ea;c > 0.035 Lq). Both quantities were derived by comparing 
the stars' location in the HR diagram to the theoretical pre-Main Sequence evolutionary 
tracks for Z — 0.3 Zq of Siess et al (1997). The luminosities and temperatures of the 
stars were computed excluding the F336W magnitude from the fit, since it is affected by 
non-photospheric emission. 
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Fig. 9. — Mass accretion rate as a function of age for Classical T Tauri stars in different 
star-forming regions (adapted from Muzerolle et al 2000). The solid line represents the 
fiducial model from the viscous disk similarity solutions of Hartman et al (1998). The 
errorbar indicates, for the Galactic stars, typical uncertainties in ages from the HR diagram 
and typical uncertainties in accretion rates due to variability. Our result for the field of 
SN 1987A in the LMC is marked with a black square. Its uncertainties are discussed in 
section 5. 
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Table 1: Log of the observations centered on Supernova 1987A. 



Filter Name 



Exposure Tiuie (secouds) 



Comments 



September 1994^ February 1996'^ 



July 1997^ 



F255W 
F336W 
F439W 
F555W 
F675W 
F814W 



2x600 2x600 

2x400 350+600 

2x300 2x300 

2x300 2x300 

2x300 2x300 



2x900 1100+1400 



2x1300 
2x800 
2x400 
2x300 
2x300 
2x400 



UV Filter 
U Filter 
B Filter 
V Filter 
R Filter 
I Filter 



F656N 



1100+1300"^ 



4x1400 



Ha 



"September 24, 1994, proposal number 5753 (SINS collaboration, PI R.P. Kirshner). 
''February 6, 1996, proposal mmibcr 6020 (SINS collaboration, PI R.P. Kirshner). 

'^July 10, 1997, except for F502N taken on July 12, 1997, proposal number 6437 (SINS collaboration, PI R.P. 
Kirshner) . 

'^February 3, 1994 proposal number 5203 (PI J. Trauger). 



